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Abstract 

Ov 

' The color dielectric modification of the gluon dispersion relation in a dense QCD medium 

suppresses both the soft and collinear gluon radiation associated with jet production. We compute 
the longitudinal and transverse plasmon contributions to the zeroth order in opacity radiative 
energy loss in the 1-loop HTL approximation. This is QCD analog of the Ter-Mikayelian effect in 
QED and leads to ~ 30% reduction of the energy loss of high transverse momentum charm quarks 
, produced in a QCD plasma with a characteristic Debye mass fj, ~ 0.5 GeV. 

m ■ 

o : 

m ■ 1 Introduction 

O ; 

Jet tomography in ultra-relativistic nuclear collisions can be used to map out the density of the 
produced QCD plasma from the suppression pattern of high transverse momentum hadrons EQ-[H]. 
The quenching of jets |10Ullj is mainly due to the medium induced radiative energy loss of high energy 
partons propagating through ultra-dense QCD matter. However, even if final state multiple elastic 
and inelastic interactions are neglected, gluon radiation associated with hard QCD processes already 
.^J | softens considerably the lowest order jet spectra. In elementary particle collisions, this associated 

radiation can be taken into account through the Q 2 (DGLAP) evolution of the hadronic fragmentation 
functions. In a plasma, even the associated radiation is modified by the dielectric properties of that 
medium as first pointed out by Ter-Mikayelian |121 113j . The non-abelian QCD analog of the Ter- 
Mikayelian plasmon effect is the subject of this paper. A summary of our main results has been 
published in paper I|14j. In this paper, the details of Ter-Mikayelian calculations are presented. 

This work is motivated by the surprising observation of PHENIX ^ that "prompt" single electron 
spectrum from open charm production in Au+Au collisions at y/s = 130 AGeV shows no sign of heavy 
quark energy loss [!£]. In contrast, a dramatic suppression (by a factor ~ 5) was observed at RHIC 
for pions with pj- > 3 GeV originating from the fragmentation of light quark and gluons |17|-|25|. 
See a recent reviews of light quark and gluon tomography in A+A in |26j-|27j. The suppression of 
light hadrons is consistent with the expected large radiative induced energy loss of light quark and 
gluon jets in an ultra-dense plasma of density approximately 100 times higher than in ground state 
nuclei 2 . 

The spectrum of induced radiation depends on the optical thickness or opacity x = L/\ of the 
plasma, and has been computed to arbitrary order in \ n f° r massless partons in GLV [3J [2Sj . It was 
expected |28|-|30| that a similar quenching pattern should also be observed for heavy quark (c or b) 
jet fragmentation. However, in |31j it was pointed out that the a large quark mass would lead to a 
"dead cone" effect, reducing induced radiation inside the cone < M/E and that this should reduce 
radiative energy loss of heavy quarks as compared to light partons. Numerical estimates indicated 
that the quenching of heavy charm quarks may be only about a half that of light quarks. PHENIX 
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data, however, suggest that the charm quark energy loss could be even smaller than that. In I ^1] we 
showed that the apparent null effect observed for heavy quark energy loss via single electrons may in 
part be due to a further reduction of the leading order 0(x ) associated. In this sense, there are two 
opposing medium effects: (1) at 0(x ) t ne Ter-Mikayelian (plasmon) effect that arises at 1-HTL loop 
level the reduces the associated energy loss and (2) at 0(x n - 1 ) the induced radiative energy loss that 
increases the total radiative energy loss. Our method of calculating the Ter-Mikayelian effect in QCD 
is discussed in detail below. 

In the absence of radiation below a plasma frequency cutoff was estimated to reduce the 
induced energy loss by only ~ 10%. However, the Ter-Mikayelian effect on the zeroth order in opacity 
(L/A)° associated radiation was not considered up to now. The first estimates of the influence of a 
plasma frequency cutoff in QCD plasmas were reported in |32[ I33| using a constant plasmon mass 
ojq [H^-jSni- The k dependence of the gluon self energies and the magnitude of longitudinal radiation 
were not investigated in that work. We extend those results by taking both longitudinal as well as 
transverse modes consistently into account via the frequency and wavenumber dependent hard thermal 
(1-loop HTL) self energy W(oj,k) ISZj-jH]. 

As noted in I the dielectric properties of an isotropic plasma lead to a transverse gluon self 
energy Hr(w,k) with 11^(^(0), 0) = 0^(0) ~ ^ 2 /3, where \x ~ gT is the Debye screening mass of a 
plasma at temperature T in the deconfined phase. In addition, long wavelength collective longitudinal 
gluon modes arise with 11^(^(0), 0) = u^iO). This dynamical gluon mass, ^/(tc&(A;) — k 2 ), suppresses 

not only the radiation of sub-plasmon co < w p ;(|k|) gluons but also shields against the collinear k± — > 
singularities. 

In this paper we study under what conditions the assumption of using a k independent effec- 
tive plasmon mass and neglecting longitudinal modes may be adequate. We show below that this 
simplifying assumption is surprisingly accurate (~ 10%) if the asymptotic mass HO], = uie/V^ 
rather than the k = ui p i = ^/2/3woo is employed. The accuracy of the approximation also improves 
dramatically as the mass of the heavy quark increases. 



2 Jet Production in the Vacuum 

In order to introduce the method that we use to compute the QCD analog of the Ter-Mikayelian |12| - 
|13j effect on the zeroth order in opacity radiation, we consider first the well known case of radiation 
in the vacuum. 



J(P) J*(P) 

FIG. 1 Illustration of a bare jet vacuum to vacuum amplitude. 

The vacuum to vacuum amplitude, illustrated in Fig.l, is given by 

- iMo = ^J ^\J(p)fA M (p) (1) 

where Am{p) = {p 2 + M 2 + ie) _1 is the jet propagator for a spinless parton of mass M in the dj 
dimensional representation of SU(N C ). We will ignore spin effects throughout since they are irrelevant 
in the soft radiation limit. The effective jet source current, J, creates an invariant distribution of jets 
as given by 

where E p d 3 Nj/d 3 p = dj\ J(E p , p)| 2 /(2(2vr) 3 ) with E 2 = M 2 + p 2 . 
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FIG. 2 The one gluon associated radiation amplitude in the vacuum. 

The one gluon radiative correction amplitude, iM\ , to the jet spectrum in the vacuum is illustrated 
on Fig. 2 

_ iMl = -^^ J | J(p + k)\ 2 A M (p + k) 2 A M (p)(2p + kyD$(k)(2p + ky. (3) 

The free gluon propagator in the axial (u^A^ = 0) gauge is 

where m is an infrared mass in the massive gluon scheme introduced by need to study the m = 
limit. In the next section, the dynamic polarization tensor will replace m. 

The imaginary part of this amplitude contains real and virtual radiation corrections to iMq. The 
inelastic one gluon radiation contribution is obtained by the Cutkosky rule [Hj 

A M (p)A m (k) -> {-2TTi) 2 5{k 2 -m 2 )5(p 2 -M 2 ) . (5) 

This gives a contribution the jet plus gluon production rate 



2fcnMi| rad = I ^5(p 2 -M 2 )djJ^5(k 2 -m 2 )\J(p + k)f 



C >(Q 2 -M 2 )^^-(Q 2 -M 2 ) 2 -4^-4M 2 + m 2 ), (6) 



(Q 2 - M 2 ) 2 [ ' (uk) K ^ ' (uk) 

where 

Q 2 = {p + kf = 2(pk) + k 2 + M 2 . (7) 

In general, the result depends on the gauge parameter because the external color source J 
breaks gauge invariance. However, in the soft gluon limit, k± <C k + <C p + , where k±_ is measured 
relative to the p axis of the jet, one can extract the familiar DGLAP soft radiation spectrum with any 

for which (uk)/(up) ~ k + /p + = x <C 1. For typical light cone kinematics of interest, the momenta 



3 



are expressed as 



k» = [xE + ,{kl+m 2 )/xE + ,k ± ] 
p» = [(l-x)E + ,(M 2 + kl)/(l-x)E + ,-k ± ] 
k\ + (1 - x) 2 m 2 + x 2 M 2 



2pk 



x(l — x) 



^2 , 2 „ ,2 k 2 + (1 - x)m 2 + xM 2 

Q 2 = pk + m 2 + M 2 = -± v - ' . 8 

x{l — X) 

Note that the vertex factor {...} in Eq. © reduces to (2pk + m 2 )/x in the x — > limit in both the 
A + = light cone and temporal A = gauges. 

We assume that J is slowly varying J(p + k) ~ </(p) for soft radiation, and therefore, in soft 
radiation approximation the spectrum can be extracted as 

2ImMi| md « | d 3 iVj | d 3 iVj°) , (9) 

leading to the finite mass generalization of the small x invariant DGLAP radiation spectrum 

dN^ _ dNP Cja s _ Cja s 1 

W d 3 k ~ X dxd 2 k ± ~ xn 2 (Q 2 -M 2 ) vr 2 k\ + (1 - x)m 2 + x 2 M 2 ' 1 j 

3 Ter-Mikayelian effect 

In this section we want to compute (zeroth order in opacity) associated radiative quark energy loss. 
In soft gluon limit, the result should not depend on the choice of gauge, as long as (uk)/(up) « x <C 1 
is satisfied. We simplify our calculations by choosing the temporal axial gauge. However, we have 
explicitly checked that, in soft gluon limit, the same result is obtained using the light cone gauge. 

The radiative heavy quark energy loss in hot dense medium involves both transverse and longitu- 
dinal gluon radiation. In temporal axial gauge gluon propagator has the following form (see appendix 
A): 

n 5^ 1 1 1 \ 



io 2 e T - k w e L 



where transverse [Pfw) and longitudinal (Q^u) projectors are given in terms of g^ v = g^ u — and 
% = g^ u Uu, by 

r iiv — y^ v —2 ' \ > 

^ _ l Ufj,k v + k^Uy iP'k^ky , {u • k) 2 . . 

^ " ^ (u-k) + Kfc) 2 " ^ j ~fc^2- W 

and 

are transverse and longitudinal dielectric functions respectively |1T] . 
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In order to compute the associated radiative energy loss when the hard process is embedded in a 
dielectric medium, we need to compute the squared amplitude of Feynman diagram diagram |M TOC ;|, 
which represents the source J that produces an off-shell jet with momentum p', which subsequently 
radiates an gluon with obeying the dispersion relation, u(k), of the medium with momentum k. The 
jet emerges with momentum p and mass M. Since our focus is on heavy quarks, we neglect the thermal 
shifts of the heavy quark mass here. To calculate this rate we use the optical theorem via 

'^'' ApffL - 2imM ™ - / <i6> 

where, in axial gauge, ItuMtm-, is the imaginary part one Hard Thermal Loop amplitude of the the 
diagram cut as shown in Fog. 3. 

ItuMtm 
I 




J(p+k) p I J*(p+k) 



FIG. 3 shows the cut amplitude that contributes to the associated radiative energy loss in the 
medium . The dashed line shows which propagators are to be put on-shell with respect to the 
dispersion in the medium. The blob denotes the resumed HTL self energy using Eq. (|lll 1141 115|) . 

We assume, as in [I], that J varies slowly with p, i.e. that J{p + k) « J{p) and neglect high x 
spin effects, amplitude in Fig.3 before cutting is 

The contribution to medium on-shell radiation is 

/cfiv d^k 1 

(2vr)32£ 1 J{P)]2 (2^* (Q 2 - {2p + k ^ D ^ 2p + ^ ^ 

where we have used jlH] l/(p 2 — M 2 + ie) — > (—2iri)5(p 2 — M 2 ) for the cut quark propagator in Fig.3 
(neglecting its small thermal shift). As shown in Appendix B„ the contour integration over lo for 
on-shell quark and glue is equivalent to replacing D^ u by: 

iV - -PA- 2 ™)^ 2 tT - k 2 ) - ^(-2ni)5(e L ), (19) 
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where we only keep positive energy contribution. A small 0(m g /E) quasi elastic contribution is 
neglected in this approximation as discussed in App. B. We obtain in this way 

d*Nj [ £N™ = [ *£-tfH|J(p)| a / C R g 2 1 



(2tt) 3 2£ ^" J s (2tt) 3 (Q 2 - M 2 ) 2 
x(2p + kr(-P^5(u 2 e T - k 2 ) - 9^S(e L ))(2p + k) v (20) 



It is convenient to choose coordinates such that p = (E, |p|, 0, 0), where E = y/p 2 + M 2 , and 
k = (uj, |k| cos 9, |k| sin 6,0). With this kinematics we get: 



(2p + k) u P up (2p + jfe)" = -(2p + fc)*(<5y - ^)(2p + £:) j = -4p 2 sin 2 

A; 2 



(2p + k) v Q up (2p + = -(2p + fc)*(^)(2p + = -(2|p| cose + |k|) 2 . (21) 
The integrated radiation yield than becomes: 

, C R a s 1 



J d 3 Nj M = J dud cos 0k 2 d\k\ 



vr (Q 2 - M 2 ) 2 



x {4p 2 sin 2 05{u 2 e T - k 2 ) + (2| P [ C ° S \ + |k|)2 5(e L )} (22) 
where Q 2 - M 2 = w 2 - k 2 + 2(Euj - |p||k| cos (9), as in Eq. (0). 

The Ter-Mikayelian modified associated (0 th order in opacity) radiative energy loss spectrum is 
defined as dl g = uod^Ng™ . Using this, the transverse and longitudinal contribution to the th order 
radiated energy loss per wave number is given by the following: 

dI T C F 4k 2 p 2 o; 2 ( W 2 -k 2 ) /-i a s (.Q 2 - M 2 ) . 2fl 

' d cos , —7: ttttt^- sm ft 



d|k| vr w 2 ^2 _ ( w 2 _ £ 2 )2 7 (Q 2 -M 2 ) 

dI L C F 4k 2 p 2 K - k 2 ) a5(Q 2 ~M 2 ) / k 2 \ 2 

ucosy — — 2 2 2 cosfl + (23) 



d|k| VT ^2 _ ((J 2 _ £ 2 ) 7 (Q2 _ M 2)2 ^ 2 £2 

where we keep only the forward, > 0, emission to isolate the energy loss of the nearside jet. In Eq. 
l(2*3*|) ujt and Wi are positive zeros of (lu 2 €t — k 2 ) and el respectively. The angular integration can be 

performed analytically if as does not run, but it is not particularly useful. We perform the integration 
using no momentum cutoffs and running coupling according to the "Frozen a model" |46j 

A/TT 

as(Q 2 -M 2 ) = Min{0.5 , q2 _ m2 J > (24) 

QCD 

where (5q = ^ for effective number of flavors rif « 2.5 and Aqcd ~ 0.2 GeV. Figures 4 and 5 show 

transverse an longitudinal contribution to the charm and bottom radiative energy loss at different 
Debye masses. 
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FIG. 4 The th order in opacity contribution to Charm quark radiated energy loss for 15 GeV jet 
is shown as a function of wave number. The dashed-dotted curve shows what would the energy 
loss be if gluons were treated as massless and transversely polarized. From top to bottom (left to 
right) solid (dashed) curves show medium modified transverse (longitudinal) contribution to the 
energy loss for Debye mass 0.25 GeV, 0.5 GeV and 1 GeV respectively. 
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FIG. 5 The th order in opacity contribution to Bottom quark radiated energy loss for 15 GeV jet 
is shown as a function of wave number. The dashed-dotted curve shows what would the energy 
loss be if gluons were treated as massless and transversely polarized. From top to bottom (left to 
right) solid (dashed) curves show medium modified transverse (longitudinal) contribution to the 
energy loss for Debye mass 0.25 GeV, 0.5 GeV and 1 GeV respectively. 

We see that longitudinal contribution to the energy loss is significant only in the u> region around 
Debye mass. As temperature T of the medium increases the transverse contribution to the energy loss 
decreases, while longitudinal contribution increases. Therefore, we can conclude that for lower Debye 
masses (lower temperatures) the longitudinal contribution to the energy loss is negligible, but at high 
enough temperatures it may become comparable to the transverse contribution. 

Fig.6 shows the integrated th order fractional energy loss for charm and bottom quarks in a hot 
plasma with Debye mass /U = 0.5 GeV. 



Charm 
M=1.5 GeV 




Bottom 
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5 10 15 20 25 
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FIG. 6 The th order in opacity fractional energy loss for charm and bottom quarks in hot medium 
with Debye mass fj, = 0.5 GeV, and zero momentum cutoff is shown as a function of the charm 
quark energy. The upper (lower) solid curve shows transverse fractional energy loss for charm 
(bottom) quark. The dashed curves show the negligible additional effect of longitudinal plasmons. 



We see that the longitudinal plasmon contribution is indeed negligible for both charm and bottom 
quarks. 

In order to study in more detail how the specific 1-loop HTL dispersion relation affects the results 
we show in Fig. 7 the dynamical gluon mass as a function of wavenumber for different screening scales. 
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FIG. 7 One loop transverse plasmon mass m g (|k|) = Vuj 2 - k 2 is shown as a function of quark 
momentum |k| 2 . We see that m g starts with the value io p i = fi/V3 at low |k|, and that as |k| 
grows, m g asymptotically approaches the value of = (J,/V% in agreement with |4l)j. 



The main feature to note is the relatively rapid transition from the static uj p \ = {i /\^3 to an asymptotic 



w-oo = as emphasized in [40] . The effect of replacing the dynamical mass by a single effective 

value, either lu p i, moo, or fi is shown in Fig.8. 
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FIG. 8 Medium modified zeroth order in opacity charm quark fractional energy loss is shown as 
a function of quark energy. Full curve shows the transverse energy loss using Eq. J23J. Dot-dot- 
dashed, dashed and dot-dashed curves show what would be the transverse energy loss if we define 
gluon mass as to p i, mi n f ty and /i respectively. 

Figs. 7,8 demonstrate that a remarkably good approximation to the the Ter-Mikayelian effect can be 
obtained by approximating m q {k) g mpp (see also Appendix Bl). 
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FIG. 9 The reduction of the fractional associated vacuum energy loss for charm and bottom quark 
due to the QCD Ter-Mikayelian plasmon effect is shown as a function of the quark energy. The 
upper (lower) dashed-dotted curve shows the vacuum energy loss for charm (bottom) quark if 
gluons are treated as massless and transversely polarized. The upper (lower) solid curve shows the 
effect of using the medium dispersion for gluons. We take fj, = GeV for the vacuum and [i = 0.5 
GeV for the medium cases. 

In Fig. 9 we compare the associated energy loss in the vacuum and a medium with fj, = 0.5 GeV. We 
see that for charm quark medium energy loss is significantly reduced by ~ 30%. The Ter-Mikayelian 
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plasmon effect therefore enhances the yield of high transverse charm quarks relative to the vacuum 
case. On the other hand, we see that for bottom quark Ter-Mikayelian effect is negligible as is the 
absolute associated radiated energy because the dead cone jSJ is so wide in that case. 
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FIG. 10 Medium modified zeroth order in opacity fractional energy loss is shown as a function 
of Debye mass /i for different mass quarks. We see that medium effect is more important for 
light (dot-dashed curve) and charm (dashed curve) quark energy loss, since it leads to the strong 
suppression of the associated radiation. On the other hand, bottom quark (dot-dot-dashed curve) 
energy loss shows only a weak dependence on the medium. 

So far, we have concentrated on the Ter-Mikayelian effect for only heavy quarks. The generalization 
of the plasmon effect to light quarks is not trivial due to the fact that the light quark vacuum energy 
loss is infrared divergent. However, as shown in Fig. 10 we see that in a QCD medium, dynamic 
polarization naturally regulates infrared divergences for light quarks, since both quarks |47| I48j and 
gluons acquire a finite self energy. Confinement in the vacuum naturally limits the effective screening 
scale to fi V ac > ^qcd- 

4 Conclusion 

In this paper we computed both the transverse and longitudinal contributions to the lowest order in 
opacity radiative quark energy loss. We have shown that longitudinal contribution can be neglected 
for the energy range of experimental interest. We have also seen that transverse polarization can be, 
for moderate range of temperatures (0.5 < fi < 1 GeV), approximated by simple form D^ v — k f^i , 

where is transverse projector and m g « moo = It is remarkable how well the effect of 

medium polarization 11^(0;, k) can be approximated numerically in this simple way. In a subsequent 
paper |44j . we will use this approximation to simplify the calculations of medium induced radiative 
energy loss. 

In the appendix B it was also seen that the poles of the gluon propagator give both the contribution 
from quasi-particles and "particle hole excitations" (which correspond to the energy u smaller than 
momentum |k|). However, we have shown that the particle hole contribution to the radiative energy 
loss is negligible compared to the radiative one. 

The next step beyond the 1 loop HTL approximation employed here will be to compute the 2 loop 
(1st order in opacity corrections) to the associated energy loss. Cutting two HTL loops diagrams, we 
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pick up three different contributions to the energy loss. Two of them give corrections to the plasmon 
mass and particle hole energy loss respectively. These are perturbative higher order corrections to 
the results here. However, a third contribution, in which the poles of two HTL propagators give one 
quasi-particle and one particle- hole excitation give rise to additional energy loss. This contribution 
corresponds to the generalization of the first order opacity correction to medium induced radiative 
energy loss computed for light partons in [I]. The results were summarized in I |14j . and the details 
of that computation will be reported in a subsequent paper 



Acknowledgments: Valuable discussions with I. Vitev, Z. Lin, J. Nagle, and W. Zajc on heavy 
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Office of Science, Office of High Energy and Nuclear Physics, Division of Nuclear Physics, of the U.S. 
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A Gluon propagator in temporal axial gauge 

We recall here some of the basic tensorial properties of the gluon propagator, in axial gauges (u tM A fJi (x, t) = 
0, u M a fixed four vector). 

The transverse projector P^ v (with respect to both k^ and u^), is given by Eq. (|12|) . Note that 
u 2 = u 2 — (ku) 2 /k 2 is well defined even in light cone (u 2 = 0) gauge. 

The orthogonal longitudinal projector, Q^ v (with respect to k^), is given by 



Ufj,ku + kpUy u k^k v 



^ " [9 ^ V (u-k) + (u-k) 2 ] ~ k 2 U 2 

/ll^ll 11 \ 1 

= (ukf + k»k v u 2 - (u»k v + Wu v ) (uk) 919 . 7 ^ (25) 

V u z ) u z k z — [uk) z 

In terms of these projectors the free gluon propagator in the axial gauge is 



D^ v (k) = -A(k)(P^ u + (3(k,u)Q^) 



1 f r _u^ + k^ +u2 k^ 



k 2 + ie V (uk) (uk) 

where 



0M = 7^2 • ( 2 7) 



k 2 u 2 
(ukf 

is a kinematic factor. In the temporal axial gauge, u = (1, 0, 0, 0), (3 = 1 — k 2 /lj 2 and the projectors 



reduce to 



Pfiv = Qvv = 0, if u, v = 

Pij = -Sij + Qij = if i,j = 1,2,3. (28) 

k k 

In a medium with four velocity u^ the gluon acquires a temperature dependent self energy Ii^ v 
in addition to its vacuum self energy. The one- loop (Hard Thermal Loop [lUj (HTL) or equivalently, 
eikonal linear response [l^) self energy can be decomposed as 
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IV = n L iv + u T p^ (29) 

where 

R/J.1/ = "=2"" (^) 

is the longitudinal projector in covariant gauges with respect to ■u' 4 but transverse with respect 
to k^ . Therefore, the HTL self energy is transverse with respect to k, kU = lik = 0. Note that 
RP = PR = 0, R 2 = R, and 

QRQ = Q/(3{k,u) . (31) 

where (3 is given by (|27|) . 

The HTL medium modified gluon propagator (D^ u ) can be obtained by solving the Dyson equa- 
tion @7] 



= D L Qr + n r P^ = D$-D$II* ) D( iv 

(32) 



uj 2 e T - k 2 u 2 e L 

Note that, in general |41| . we would have one extra term r]^^- in the gluon propagator, where rj 
is a gauge parameter. However, in temporal axial gauge |47j . 
= 0, if fi, v = 0. 

Therefore, in temporal axial gauge r\ has to be equal to zero, i.e. extra term vanishes. 

B Cutkosky rules for gluon propagator in the medium 

In this appendix we want to justify using the Cutkosky rules for gluon propagator in the medium, 
which is represented by Eq. (|19j) in our computations. 



Diagram M 




J(pj+k) p, J*(pj+k) 

ro FIG. 11 

corresponds to M n , where M n is the amplitude of the following diagram: 
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{iM n )=Tr{^j — — \ J{p +k)\ ((pj + ^2 _ M 2 + ie) 2 p 2 _ M 2 +ie 

x(-ig s (2pj + fc)^r ew )^^i7r^w^^i7r^...i7r^3^(-i^(2p J + fc)^T ao )X33) 

/fr + Ar + fe z + le 

where w^Qc) is the amplitude of the diagram shown on Fig. 13, averaged over thermal momentum 
distribution n eq (p) given by Eq. (|35|). 



iTr^(fc) = J d 4 pn eq (p){ - - - - (-ig s (2p + ky){-ig s {2p + kf) + 

+7 irj-rHl/.(2p " fe) A ')(-»5.(2p - AO") + 2^}Ei T « r W' t ( 34 ) 

where ^ - jT a T\,j^ = ^5^, and n eq (p) is the equilibrium momentum distribution 41 at tempera- 
ture T including both quarks and gluons 
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1 

2 

Here 



n eq (p) = ^(Qt q + Qe q ) + NG eq . (35) 



INt rP 

Q% = j^29(±p°)5(p 2 )(eM±Y) + ir ' (36) 

and 

2N r>° 

G eg = j^29(p°)5(p 2 )(eMY) " 1} (37) 

are quark (antiquark) and gluon distributions respectively. Nf is the number of flavors, and N is 
the number of colors. 

In the soft gluon limit tt^ v becomes 

* =- 9 °J d P MTTe UM (38) 

in agreement with 41 . 

It is easy to prove that Tr^ u is transverse, i.e. TT^ v k y = k^ir^ = 0. 
Therefore, in covariant gauge ir^ can be decomposed: 



= U T P^ U + U L R^, (39) 
where P^ u and are given by Eqs. (|12[) and ()3Uj) respectively. 
Using Eqs. (3H1 EH dl USB we get 

U L = = - 9 l ^/ n ^ ^ p ^) WW-(P")(^) (40) 



u 2 A; 2 — (ku) 2 J (pk) 



U T - ^ _ / ^ { (^(-JJM^ _ „*) + o*)^) _ »g>} (41) 



By replacing n eq {jp) from Eq. (|33|) we get 



(lo 2 - k 2 )// 2 . x , ^' - jk 



and 



k 2 2 k w+k 



nT = | + Mz±V (1+ « log |^|) (43) 
2 2k 2 2 k cj+ k 



where /i = g s T\jl + . 

These results are in agreement with |41j . and they lead to and ej- which we use in this paper. 



14 



Using 7r Atiy from Eq. M finally becomes: 



m = H)/^MI^ p2j _^ + ie 



x / C R gl 



d 4 k 



1 



where D 



ft]' 



(2vr) 4 (Oj + A;) 2 - M 2 + ie) 2 
fc^_n T +ie ~~ an d we assume that e is positive. 



(2pj + kYD tlu {2p J + k) u 



(44) 



Lets now compute 



((PJ + k) 2 - M 2 + ie) 2 



(2 PJ + kfD^pj + k) v = J M 



fripj,k) fi{pj,k) 



k 2 -U T + ie k 2 -Il L + ie 



)] (45) 



Since we are interested only in the radiative energy loss, we assume that initial jet is off-shell, and 
therefore we have that /t(pj, k) and fi(pj, k) are analytic functions. 

^From Eqs. (jiUl l4*T|) we see that Hr and Ul can be written as J2 P ^^(pk)^ where £,t(L)(Pi k) are 
analytic functions, and where we have assumed that we have discrete set of p's. 

Since discussion for both parts of the integral ()45j) is the same, we will consider only one part: 

f(pj,k) 



dtu- 



u 2 - k 2 - £ 



g(p,fc) 

P (pk) 



+ ie 



(46) 



Note that ^j^y is infinite when u = 



For simplicity, suppose that all ^ are different. We can order them in such a way that uj\ < L02 < 
< uj n < |k|. 

Than, for fixed |k| zeros, ±(0^— ie), of (uj 2 — k 2 — V p ^fe^ + ie) can be found graphically: 

co 2 -k 2 




FIG. 14 shows schematically the zeros Qi for lo 2 — k 2 — ^2 p ^fcffi at fixed |k 
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We see that we have n' < n different solutions fij < |k|, and exactly one solution > |k|. Then, 
if we close the contour in Eq. Q46|) in clockwise direction, we will pick up only positive poles, i.e. 
(fii - ie). 



Therefore, 



/ duJ 2 & ^ Lk) ■ = ( " 27ri) E Res ^ + (-^)Res(Q ). (47) 

In the following subsection, we will test the relative magnitude of the two contributions for the 
case of the dominant transverse excitations. Solutions Oj <|k| correspond to particle hole excitation, 

and we will prove that this contribution is negligible. 

B.l Simplifying the gluon propagator in hot dense medium 

The transverse response 



obeys the sum rule |4*8] 



duujp T {uJ, |k|) = 1 (49) 



We can write 



pr(u, |k|) = 0(u, |k|) + 5{u 2 e T - k 2 ), (50) 



where first part in this formula corresponds to the particle hole excitation, and the second part 
corresponds to the delta function contribution. 
It is easy to show that 

5(u 2 e T - k 2 ) = 5(io 2 - u T )Z T {k) (51) 
where lot is the transverse plasmon spectrum, and 

2uj t (uj t -k 2 ) 

T( } ~ uJ 2 fi 2 -(uJ 2 -k 2 ) 2 ' (52) 

The sum rule reduces to 

r°° -» 

Z T (k)+ duuP(u,\k\) = l (53) 
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FIG. 15 shows Zrfjt) as a function of |k|. Full curve shows exact Zt(ojt, |k|), while dot-dashed 
curve shows Zxiyk 2 + m 2 , |k|), in which ujt is approximated by v/k 2 + m 2 . 



In Fig. 15 we see that Zt(ujt, |k|) is approximately equal to 1 for whole region of |k|. Therefore, 
we can conclude that the contribution form the particle hole excitation is negligible. 

In order to justify this conclusion we will estimate the particle hole contribution to the energy loss. 
For simplicity, we will consider quarks to be massless. 

^,From Eq. (|22j) we see that 

= / l^ricosflk 2 d|k| ggg {Q 2^ M 2 r 4p 2 sin 2 fl/^,|k|) < 

2 C R a^ X /V 2 „r, , nr ,i*iw A+ X ^ 
- - — 

where x = 

|k| 



< - - 5 / k^lkl xdx /3(x,\k\) (Ln(^— ) - 2) (54) 

E 7T J I — X 



(Ln(i±^) - 2) < (Ln(-^-) - 2) - Ln(l - x) Q(x - 0.95) (55) 
1 — x 0.05 

Using Eq. (|55|) we get: 

^E~ < 1^?- { (Ln( 0^5 } " 1} / /k 2 ^^|k|)- 

- f 1 xdxLn(l-x) [ lPl d\k\ k 2 |k|) } < 

iO.95 i0 

< |«eC { 2 /" rf |k|"- Z ^» - /' rf* Ln(l - ,) /" dtf i? |k|) }(56) 
E IT Jo 2 Jo.95 JO 

where in the first integral we have used Eq. (|53|) . 
Using the fact that |1 — Z<rW)\ ~ 9 , N we get 

2(k^+mg) 
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.max 



2<W r^^^^y^C^ 



E ir Jo vr E 2 

In the x — ► 1 region (3(x, |k|) |48| can be approximated by 

m 2 (1 — x) 

(3(x, |k|) « ^ 9 V i (58) 

(2k 2 (l-x)+m 2 ) 2 

Then, the second integral in Eq. Q56JI can be performed analytically, leading to: 

We see that the particle hole contribution to the fractional energy loss is on the order (m g /E). 
Since m g <C E we can conclude that this contribution is negligible. 



Finally, on Fig. 7 we saw that lot ~ y k 2 + m 2 . Dot-dashed curve on Fig. 15 represents Zt("J k 2 + m 2 , 
as a function of |k|. Again, with the exception of the small |k| < 1 GeV region, we see that 



.^(y'k 2 + m 2 , |k|) is approximately equal to 1 for whole region of |k|. Thus, we can approximate 

<5(cA T - k 2 ) « 5{oo 2 - (k 2 + m 2 )). (60) 

Since longitudinal contribution is negligible, we can conclude that a gluon propagator in a hot 
dense medium can be approximated by 

D ^ « - ,2 - ■ ( 61 ) 
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